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Using Ecological Niche Modeling for Biodiversity Conservation Guidance in the Western Podillya 
(Ukraine): Amphibians. Tytar, V., Mezhzherin, S., Sobolenko, L. — Maximum entropy niche modeling 
was employed as a tool to assess potential habitat suitability for 13 amphibian species and to map their 
potential distribution in the Western Podillya (Ukraine). The predictor variables used were of climate, 
topography and human impact (assessed by the Human Footprint). The “mean temperature of coldest 
quarter” and “isothermality” were two of the most important factors in predicting habitat suitability and 
distribution. Another profound contribution has been displayed by the Human Footprint, meaning that 
human infrastructure may benefit amphibians, a phenomenon that perhaps is much more widespread 
than thought. Areas have been distinguished that in the first place should be of interest to nature 
conservationists targeting amphibians (exemplified by Bombina variegata) and a map summarizing 
species richness was produced. 
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Mogeauposanne skOOrM4ecKOM HMMM KAK MHCIpyMeHT [JIA MIAHNMPOBAHMA Meponpuatuit, 
HalipaBeHHbIX Ha coxpamenne 6uopasHoo6pa3ua Samaquoro Ilogonba (YKpanna): ampuéun. 
Tutap B., Meaoxepun C., Co6onenxo JI. — Mogenupopanue skonormuecKoi HMM MeTOOM 
MaKCMMasIbHOM JHTPONMM ŐÞIIO MCHO7b30BaHO JIA OMEHKM ycnoBui npebbiBanna 13 BugoB ampubuit 
u KapTMpoBaHMA nx pacnpocrpanemna Ha TepputTopuu 3anagHoro Ilogomba (YKpamua). Cpem 
IpegMKTOpOB ObIIM MCHOMb3OBaHbI NOKAZATEIM KAMMaTa, pembecdba n aHTpOMOreHHOrTO BOsAelcTBMA 
(oWeHMBaeTCA NO MHTeETpMpOBaHHOMY MHJeKcy «semoBeyecKkMi creg»). Cpegu Baxneňmnx þakTOpOB, 
KOTOpble ONpegenIoT NpurogHocTh M pacnpegenenne MecT obuTaHMA aMnónÄŇ, Opin «CpegHaAA 
TeMIepaTypa caMOTO XOOJHOTO KBApTaNaA» M «M3OTEPMMYHOCTD». J|pyroli CynyecTBeHHDIŇ BKJIAJ BHOCHT 
MHJeKC «YeNOBeYecCKOTO Cega». ITO MOXKET O3HAYATD, YTO MHpactpykTypa, Co3qaHHaA YeNOBEKOM, 
copmupyer 6NAronpnATHDIe YCAOBNA PIA aMpuOMM N, BOSMOXKHO, ITO ABNEHNE GONee paclIpoctpaHeHo, 
4eM IIpefcTaBslaocb paHee. OTMEYEHDI TeppUTOpUM, KOTOpHIe B NEpBYyIO O4epe|{b JOMDKHBI MpeACTaBAATS 
wHTepec pia oxpanbi am@ubuit (Bug Bombina variegata B3aT B KayecTBe npumepa) n cosgaHa KapTa, 
KOTOpaa o606NAeT BUTOBOe GoraTcTBO MccieqOBaHHOLO pernoHa. 


Karyesse cxoBa: Maxent, mopennposanne 9KOMOrMYeCKOM HMM, MOJENN pacnpocTpanenna 
BUOB, ampuoun, Ykpanna. 


Introduction 


Amphibians are among the most threatened taxonomic groups worldwide. Numerous studies have 
documented declines in amphibian species abundance across the globe. Habitat fragmentation, degradation 
and loss, together with climate change are probably the most important drivers of population decline (Billeter et 
al., 2008). A report on the status of amphibians globally (Stuart et al., 2004) stated that about 32 % of amphibians 
are clearly threatened with extinction of which 22.5 % are poorly studied. The report also noted that over 
100 amphibians are thought to have become extinct in very recent decades and that about 43 % of all described 
species are currently experiencing population declines. Therefore amphibians represent an exceptional group 
of species that are highly sensitive to both habitat and climate change, including other factors impacting the 
environment (Beebee, Griffiths, 2005). 
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To protect amphibian species, we need a better understanding of what constitutes suitable habitat and 
where such habitats exist. Habitat suitability mapping can identify areas in need of restoration or preservation 
(Gibson et al., 2004) and guide conservation plans (Gaston, Williams, 1996). However, data on species’ distri- 
butions are often sparse, so one option to cope with this problem is to use habitat modeling approaches (Tytar, 
2011). These models describe the environmental requirements of species and use it to produce distribution 
maps that are a pivotal stage in targeting conservation and recovery efforts (Elith et al., 2006; Peterson, 2006). 
Modeling can be used for revealing species ecological requirements and relationships between the distribu- 
tion of species and predictive variables, as well as the importance of each variable in model building (Araújo, 
Guisan, 2006). 

A variety of distribution modeling methods are now available for predicting the potential geographi- 
cal range of a species. Unfortunately, the performances of most species distribution modeling methods are 
poor when sample size is small (for instance, < 10). Under these circumstances maximum entropy distribution 
(Maxent) modeling may be a good choice. Maxent, unlike other distributional modeling techniques, uses only 
presence and background data instead of presence and absence data. This method has been shown to perform 
well in comparison with alternative approaches (Elith et al., 2006). 

In this study, we employed maximum entropy niche modeling as a tool to assess potential habitat suitability 
for amphibians and to map the potential distribution of this group for the area of the Western Podillya in Ukraine. 
More specifically, our objectives were to (1) identify the factors associated with (species) habitat distribution; (2) 
predict potential distributions of the species using known presence observations; and ultimately (as the primary 
goal) 3) produce a regional map of amphibian species richness for guidance of conservation measures. 


Study region 

The Western Podillya is part of the vast Eastern European Plain, bordered by the Dniester River and 
the Carpathians in the southwest. More specifically, the region of our interest is confined to a bounding box: 
Xmin = 24.73326, Xmax = 26.79984, Ymin = 48.43340, Ymax = 50.04167 (figures in decimal degrees) and cov- 
ers an area of about 27 thousand sq. km. The average altitude is 320-350 m. The climate is Atlantic-continental 
with the annual mean temperature varying around 7.54 °C and receiving annual precipitation of the aver- 
age of 650 mm. Western Podillya is occupied by the West Forest-Steppe Zone of Ukraine and belongs to the 
Continental Biogeographical Region as defined by the European Commission and the Council of Europe for 
evaluation and assessment of nature conservation (European Environmental Agency, 2002). According to the 
global land cover database GLC2000 (Bartholome, Belward, 2005), above 84 % of the area can be classified un- 
der the category “cultivated and managed areas”, leading to continuously increasing fragmentation of habitats. 
Only little remains of natural forests stands: “tree cover” (including “broadleaved, deciduous, closed”, “needle- 
leaved, evergreen” and “mixed leaf type”) occupy just 6.2 % of the study region. However, despite the seemingly 
large human pressure on the ecosystems, the region yet retains a certain amount of “wilderness”, as evidenced 
by the Human Footprint (HF). HF has been produced through an overlay of a number of global data layers 
that represent the location of various factors presumed to exert an influence on ecosystems: human population 
distribution, urban areas, roads, navigable rivers, and various agricultural land uses. The combined influence of 
these factors yields the Human Influence Index (HII) (Sanderson et al., 2002). The HII, in turn, is normalized 
by global biomes to create the HF dataset, having values ranging from 1 to 100. For the study region the average 
HF is 40.8, wherein 71 % of the area has values below 40 and 7.2 % — below 20, meaning there may be some 
good chances for successful nature conservation ventures in the region. 


Methods 


Occurrence Data Collection and Processing 

We digitized presence survey data from Sobolenko (2010) to generate the occurrence data used in the 
modeling. Georeferencing (in OziExplorer v. 3.95.4 m) was accomplished for 413 point data obtained for 
13 species: Bombina bombina (70), B. variegata (11), Bufo bufo (19), B. viridis (18), Hyla arborea (71), Pelo- 
bates fuscus (16), Pelophylax esculentus (15), P. lessonae (20), P. ridibunda (70), Rana arvalis (23), R. tempo- 
raria (29), Lissotriton vulgaris (29) and Triturus cristatus (22). Sources for taxonomy are from the IUCN Red 
List of Threatened Species, version 2013.2 (http://www.iucnredlist.org). 


Environmental Data Collection and Processing 

In most cases environmental predictors are selected based on the availability and experience that the 
variables show correlation with the species distribution (Guisan, Zimmerman, 2000). Biotic factors, which are 
challenging to model explicitly, may nonetheless be implicitly represented in the model because they strongly 
correlate with abiotic factors (Soberón, Nakamura, 2009; Tytar, 2011). In such circumstances it is reasonable to 
assume that biotic processes that lead to the species realized distribution may be captured by the relationship 
between the environmental predictor variables of abiotic character and the modeled species’ occurrence 
patterns and it is reasonable to consider modeling the distribution only with selected environmental variables 
and meaningful climatic factors identified to be of most importance to amphibians (Girardello et al., 2009) and 
topographic features considered important for shaping amphibian communties (Ford et al., 2002). In this work, 
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we used climatic predictor data, sourced from the Worldclim dataset (Hijmans et al., 2005). The 19 Worldclim 
variables represent annual trends (e. g. mean annual temperature, annual precipitation) and extreme limiting 
environmental factors (e.g. temperature of the coldest and warmest months, precipitation of the wettest or 
driest quarter) and are known to influence species distributions. The topography predictor variables selected 
were: elevation, topographical wetness index (TWI) and aspect. These predictor variables were derived from a 
digital elevation model (DEM) distributed from the SRTM 90m Database (Jarvis et al., 2008). Using SAGA-GIS 
software (http://www.saga-gis.org), 2 layers were calculated from the DEM: aspect (a proxy for the amount of 
solar radiation on the ground surface) and the topographical wetness index (TWI), because of the important 
role played by moisture in habitat selection by amphibians (Wyman 1988). The TWI combines a measure of 
the upslope area and slope to predict the hydrology of a given location (Sorenson et al., 2005). Small values 
represent upper catenary positions (dry), and high values represent lower catenary (wet) positions. Finally, as 
a measure of anthropogenic impact, the Human Footprint (HF) data set, already mentioned above, has been 
included to the suite of variables for creating the models. All environmental data layers were spatial resolution 
rasters (~1 km) with the same extent and cell alignment, as required by most modeling software. 

Due to the high levels of correlations between many environmental variables, we filtered the initial vari- 
able set of 23 predictors based on the results of multi-colinearity analysis. High correlations among variables 
may result in highly unstable performance of the Least Squares Estimator (LSE), which will lead to problems 
for running species distribution models. Multi- colinearity can be detected by calculating the Variance Inflation 
Factor (VIF): VIF = 1/(1-R?’), where R? is the coefficient of determination. The rule of the thumb is that VIF > 10 
means multi-colinearity may influence the LSE. Multi-colinearity was analyzed by using the Multiple Regres- 
sion Tool in Statistica v.8.0. Removal of the variable with the highest values from the variable-list was followed 
by re-running colinearity diagnostics, till all the remaining values are below 10 (table 1). 


Model building and evaluation 

The Maxent software (version 3.3.3e) was utilized for modeling (http://www.cs.princeton.edu/schapire/ 
maxent/), using the default settings. Logistic output format was used to describe the probability of presence 
(Phillips and Dudik 2008), which is a continuous habitat suitability range between 0 (unsuitable) and 1 (the 
most suitable). Maxent was run ten times (using all predictor variables) for each species in order to get average 
prediction. A bootstrapping replication technique was applied to the dataset which uses all occurrence data to 
build the model. This method is optimal for dataset with few occurrences such as, for instance, B. variegata. The 
outputs in ASCII format were processed and visualized using DIVA GIS 7.5. The Jackknife analysis was used 
to indicate the most informative variables (in corresponding percentages > 10 %) and a look at the response 
curves from Maxent helps to establish the relative importance of each variable. The accuracy and performance 
of species distribution models were evaluated using threshold-independent receiver operation characteristic 
(ROC) analysis (Elith et al., 2006; Phillips et al., 2006). The area under the ROC curve (AUC) ranges between 0 
and 1. Models with an AUC value higher than 0.7 are considered acceptable (Swets, 1988). 

The logistic probabilities provide a relative indication of the likelihood of occurrence by the species, but 
they do not define predicted occurrence in the binary, presence/absence manner typically required by managers. 
Therefore, we applied three thresholds to the logistic output of each model to produce a four-category model, 
ranging from “Very Low” to “High” predicted probability of occurrence. The “Very Low” category contained 
logistical values ranging between 0 and the “Minimum Training Presence” (i. e., the logistic prediction for the 
training presence point with the lowest logistic prediction value). The “Low” category represented logistical 
values ranging from the “Minimum Training Presence” value to the “Maximum Training Sensitivity Plus 
Specificity” threshold (i. e., that threshold which maximizes the sum of sensitivity and specificity for the training 


Table 1. Remaining predictors after multi-colinearity analysis 


Taoan ga. 1. IIpequxroppi, OcTarolMeca NOCHE IpOBeqeHuA My/IbTMKOIMHeapHOroO AHAMIA 











Predictor VIF 
Bio 2: Mean Diurnal Range (Mean (period max-min)) (°C x 10) 8.3 
Bio 3: Isothermality (Bio 2/Bio 7') 6.0 
Bio 8: Mean Temperature of Wettest Quarter (°C x 10) 4.6 
Bio 11: Mean Temperature of Coldest Quarter (°C x 10) 7.0 
Bio 14: Precipitation of Driest Period (mm) 4.0 
Bio 15: Precipitation Seasonality (Coefficient of Variation) 4.4 
Bio 19: Precipitation of Coldest Quarter (mm) 4.9 
Aspect 1.0 
Elevation 5.4 
Human Footprint (HF) 1.1 
Topographical wetness index (TWI) 1.2 


‘Bio 7: Temperature Annual Range (°C x 10) 
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data). The “Moderate” category contained values ranging from the “Maximum Training Sensitivity Plus 
Specificity” threshold to the “50th Percentile Training Presence” (i. e., the threshold representing the median 
logistic prediction value for all training presences). Finally, the “High” category contained values ranging from 
the “50th Percentile Training Presence” value to 1. Binary predictions were considered choosing the “minimum 
training presence” threshold (i. e., the logistic prediction for the training presence point with the lowest logistic 
prediction value) (Liu et al., 2005). Final versions of maps were considered to benefit by a simple smoothing 
filter as one of the means for coping with observation bias (Home et al., 2007) and generalizing raster outputs: 
3 x 3 neighborhood filtering implemented in DIVA GIS was applied for this purpose. 


Results and discussion 


Factors associated with (species) habitat distribution 

The Jackknife test of variable importance showed that Bio 11 (the mean temperature 
of coldest quarter) was one of the most important factors in habitat distribution prediction 
for amphibians in the region (table 2): 12 species out of 13, particularly B. variegata (fig. 1), 
have ecological requirements dependent on Bio 11, obviously having a negative impact on 
winter survival. The only exception is P. esculentus: the modeling suggests precipitation of 
the coldest quarter (Bio19) has greater significance for the species — lower levels (around 
90-95 mm) favor predicted presence probability, whereas levels around 125-130 mm force 
the probability down to almost zero. 

The other bioclimatic variable widely affecting habitat distribution prediction for 
amphibians in the region is Bio 3 (isothermality). For at least 6 species (B. bombina, 
H. arborea, P. fuscus, P. lessonae, L. vulgaris and T. cristatus; the percent contribution of Bio 
3 > 20 %) optimality (according to the response curves) is reached at values around 26.5-27.5 
and below these there is an exclusively sharp drop of predicted presence probability towards 
zero, clearly indicating the negative influence of “temperature unevenness” over the course of 
a year (fig. 2; P. fuscus taken for an example). On the other hand, excessive isothermality too 
negatively influences the predicted presence probability, but the effect is somewhat smoother 
(in the case of P. fuscus the predicted presence probability drops to about 0.2). 


Table 2. The percent contribution of environmental variables (factors) in predicting the species 
geographic distribution models 


Ta6anya 2. Iponentupiit Bkiay pasIMTHbIX NApaMeTpoB oKpy>Katoljeit cpegbi (bakTOpos), McHomb- 
30BaHHDIX PIA MOCTpOeHMA NpPOTHOCTHYECKMX MOJelel pactipoctpaHeHMA BUOB 






































Factors associated with (species) habitat distribution > 10 % 

e) es | e@ |) se] a | S| 8 | 2] ae z 
Bombina bombina - 23:9 - 11.1 - - 16.9 - 10.1 13.5 - 
B. variegata 26.2 11.4 7 30.0 = - z = = - - 
Bufo bufo - - - 16.6 - - - 11.2 196 157 - 
B. viridis - - - 18.3 - - - 14 - 27.2 - 
Hyla arborea - 21.5 - 10.1 - - 15.7 - 10.2 14.2 - 
Pelobates fuscus 12.4 25.0 - 14.3 - - - - - 204 - 
Pelophylax esculentus - 14.0 - - - - 38.4 - - 22.4 - 
P. lessonae - 23.4 - 16.6 - - 13.2 - - 15.5 17.0 
P. ridibunda - 18.7 - 16.0 - - 14.4 - - 12.8 10.2 
Rana arvalis - 15.0 - 14.1 - - - - 18.3 12.4 - 
R. temporaria = 13.0 = 12.6 10.1 = = = 28.2 = = 
Lissotriton vulgaris - 21.4 = 12.9 = - - - - 12.9 19.4 
Triturus cristatus - 24 - 14.2 - - 15.6 - - 26.3 - 


Importance for n species: 2 11 0 12 1 0 6 2 5 11 3 
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Bio 11 


Fig. 1. Response of Bombina variegata to Bio 11: x-axis — mean temperature of coldest quarter (°C x 10); 
y- axis— logistic output (probability of presence). 


Puc. 1. Peaxyna Bombina variegata na BosgeñcrBne Bio 11: ocb x — cpegnaa TeMMepatypa HamOomee Xoro- 
HOTO KBapTazia (°C x 10); ocb y — moructuyeckum dopMata 3HaveHMi MOJENN (BepoATHOCTb IpMcyTCTBMA). 
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Fig. 2. Response of Pelobates fuscus to Bio 3: x-axis — isothermality; y-axis — logistic output (probability of 
presence). 


Puc. 2. Peaxuma Pelobates fuscus na Bospeñcrene Bio 3: ocb x — M30TEPMMUHOCTD; OCb y — NOTNCTHYeECKNĂ 
dopmata 3HadeHuit MOJENN (BepoATHOCTb IpucyTCTBMA). 


Of the non-bioclimatic variables a profound contribution in predicting the species 
geographic distribution models for 11 species has been displayed by the Human Footprint. 
The HF percent contribution reaches values ranging from 12.4 % (for R. arvalis) to 27.2 % 
(for B. viridis). Surprisingly (or not), the HF positively affects predicted presence probability 
(fig. 3; T. cristatus taken for an example). This is a highly interesting finding. Indeed, the 
relationships between human factors and biodiversity are important to assess the risk of 
extinction as human pressures are often related to large changes in biological diversity. 
However, the literature shows contradictory results. Previous studies report that human 
influence may affect species’ spatial distribution both negatively and positively (Young et 
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Fig. 3. Response of Triturus cristatus to the Human Footprint: x-axis — Human Footprint; y-axis — logistic 
output (probability of presence). 


Puc. 3. Peaxyna Triturus cristatus Ha BO3geŭcTBne þaKTOpa «4eMOBE4eCKOTO cega»: OCb X — MHJEKC «4eMOBe- 
yecKoro Ciesla»; OCb y — NOTUCTUYeCKUÄ PopMata 3Ha4veHUi MOJENN (BepOATHOCTb IPMCyTCTBMA). 


al., 2005). On one hand, human factors, such as human activities and, in particular, the 
alteration of habitats (Kiesecker et al., 2001) are major causes of biodiversity loss (Brooks 
et al., 2002). On the other hand, several studies have shown a positive relationship between 
human density and biodiversity, indicating that species-rich areas and human enterprises 
quite often co-occur (Luck, 2007). One reason may be that though human population is 
concentrated in regions critical for amphibians, there is still a substantial amount of intact 
habitat in many of these regions. Nevertheless, amphibians have been found breeding in 
a variety of habitats that are substantially different from their former pristine breeding 
habitats (Rubbo, Kiesecker, 2005), so native wildlife can often adapt to novel and altered 
habitats, given suitable conditions. In North America and in Australia, for instance, human 
infrastructure provided beneficial environments to some amphibian species (Tyler et al., 
2007). Our assumption too is that human-constructed habitats such as ponds, fish farming 
facilities etc. have realized (or are on their way to realizing?) their potential to provide 
habitat for most of the amphibian species in the study area, a phenomenon that perhaps 
is much more widespread than thought. Beyond this general trend are only 2 species: 
B. variegata and R. temporaria, but this is just because the HF percent contribution to 
the predicted presence probability of these species is not so profound — 2.4 % and 7.8 %, 
respectively. 


Potential distributions of the species 

Based on the maximum entropy modeling algorithm and using 11 environmental 
variables (VIF < 10), we obtained 13 raster outputs modeling the distribution of the 
corresponding study species. Models providing an excellent prediction have an AUC > 0.9 
(these are models for 7 species, see table 3) and fair models having an AUC between 0.7 and 
0.9 (in fact, > 0.8) have been produced for the remaining species. Models with AUC < 0.7 
are considered poor (Swets, 1988). 

The resulting models contain four categories indicating the relative likelihood of 
occurrence for each species. These categories may be used to determine whether site- 
specific surveys are needed if a management action (e. g., establishing a protected area) is 
being planned. Areas categorized as “Very Low” are the most unlikely to host populations 
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Table 3.Summary statistics for Maxent habitat suitability models 


Ta6anua 3. MTOTOBbIe cratuctuku MOpenel MpUTOMHOCTH MecTrooOuTaHNi, MOCTpOeHHbIX B IporpaM- 
Mme Maxent 














. Suitable Habitat (% of study area) 
kii Very Low Low er High oe 
Bombina bombina 72:5 9.7 11.4 6.4 0.873 + 0.015 
B. variegata 26.9 53.8 11.2 8.1 0.902 + 0.019 
Bufo bufo 61.8 22.8 5.5 9.9 0.898 + 0.036 
B. viridis 44.3 43.0 1.7 11.0 0.896 + 0.021 
Hyla arborea 17.9 61.8 11.4 9.0 0.866 + 0.022 
Pelobates fuscus 68.5 17.7 7.3 6.5 0.941 + 0.016 
Pelophylax esculentus 55.6 26.6 6.9 10.9 0.941 + 0.015 
P. lessonae 67.2 23.5 5.4 39 0.941 + 0.017 
P. ridibunda 31.3 51.7 9.5 7.6 0.870 + 0.020 
Rana arvalis 36.1 54.6 2.3 6.9 0.884 + 0.030 
R. temporaria 45.4 45.5 3.5 5.6 0.913 + 0.023 
Lissotriton vulgaris 45.7 42.6 7.3 4.3 0.917 + 0.017 
Triturus cristatus 71.4 13.5 10:9 4.1 0.937 + 0.020 


SD! — standard deviation 





kilometers 








Fig. 4. Two upper categories (“Moderate” and “High”) collapsed to identify areas of predicted presence (dark 
gray shading) for B. variegata in the study area. 


Puc. 4. [Be sepxune Kareropun («YMepeHHblii» n «Bbicoknň») OO beAMHEHDI B eMAX BÞIABIEHNA OOMacTeM 
IIpOrHo3sMpyeMoro IIpMcyTcTBMA (TEMHO-CepbIit OTTEHOK) pia B. variegata B nccnepyemom pañone. 
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of the species, and may suggest that site surveys for the species are not warranted within an 
area prior to management activities. Conversely, areas mapped as “Moderate” or “High” 
are likely very suitable for the species and suggest that surveys should be conducted prior 
to management actions to determine whether the species is present and the degree to which 
it may be impacteThe case of B. variegata can exemplify this approach. This toad has been 
considered as far as the northern margin of the home range of the species runs through 
the study area, meaning populations here are highly fragmented and more vulnerable to 
impact than elsewhere (Sobolenko, 2010). Results for B. variegata predict that only 8.1 % of 
the area is of “High” and 11.2 % of “Moderate” suitability for the species (table 3). Together 
these areas (likely to be very suitable for the species) are shaded in dark gray on the map 
(fig. 4.), and they in the first place should be of interest to nature conservationists targeting 
the species. For instance, in Ternopil Oblast the most promising areas for protection of the 
species are located alongside the Dniester River, particularly within the districts (rayons) of 
Zalischyky and Borshchiv. 

By overlaying the binary maps derived for separate species that indicate either 
presence or absence, a summarizing species richness map was produced (fig. 5): light gray 
areas are predicted to accommodate 1 to 5 species, darker ones — from 6 to 10 species, 
and the darkest — 11 to 13 amphibian species. As seen, some of the richest areas in term 
of amphibian species composition are rayons of Chortkiv, Zalischyky, Borshchiv and 
Kamyanets-Podilskiy. 

Together, in such a way predictive distribution models can be used to protect rare 
species and species’ assemblages, but as far as amphibians represent an exceptional group 
of species that are highly sensitive to environmental change it may be that protecting 
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Fig. 5. Summarized species richness map (see text for explanation). 


Puc. 5. Kapta, o6o6marmas BUOBOe HoraTcTBO MCCIEȚOBAHHOTO perNOHA (CM. TeKCT IA OObACHEHMA). 
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amphibians the human society can succeed in protecting itself from habitat degradation 
and the collapse of supporting ecosystems. 
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